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Abstract 

Inelastic neutron scattering was used to study the wavevector- and 
frequency-dependent magnetic fluctuations in single crystals of superconduct- 
ing YBa2Cu306+a;. The spectra contain several important features, including 
a gap in the superconducting state, a pseudogap in the normal state and 
the much-discussed resonance peak. The appearance of the pseudogap de- 
termined from transport and nuclear resonance coincides with formation of 
the resonance in the magnetic excitations. The exchange energy associated 
with the resonance has the temperature and doping dependences as well as 
the magnitude to describe approximately the electronic specific heat near the 
superconducting transition temperature (T c ). 



The parent compounds of the high-transition-temperature (high-T c ) cuprate supercon- 
ductors are antiferromagnetic (AF) insulators characterised by a simple doubling of the 
crystallographic unit cells in the CuC>2 planes. When holes are doped into these planes, 
the long-range AF ordered phase disappears and the lamellar copper oxide materials be- 
come metallic and superconducting with persistent short-range AF spin correlations (fluc- 
tuations). Although spin fluctuations in cuprate superconductors are observed for materials 
such as YBa2Cu306+2 [denote as (123)06+J at all hole doping levels, x the role of such 
fluctuations in the pairing mechanism and superconductivity is still a subject of controversy 
0. The most prominent feature in the fluctuations for (123)06+ x is a sharp resonance which 
for highly doped compositions appears below the superconducting transition temperature 
T c at an energy of 41 meV |l|-|3|. We show that the temperature-dependent resonance inten- 
sity is correlated with the electronic part of the specific heat of (123)0 6+:c , C el (x,T) 
including the extended fluctuation regime ("pseudogap phase"), and that the pseudogap 
temperature T* determined from the onset of the resonance intensity agrees with results of 
transport and nuclear magnetic resonance (NMR) techniques. By making absolute intensity 
measurements of the magnetic resonance, S res (q,uj), we estimate the contribution of its ex- 
change energy [0,0 to the C el (x, T) anomaly around T c . We find that the temperature and 
doping dependence of the resonance exchange energy can account for the C el (x, T) anomaly 
H, suggesting that a large part of this anomaly is due to spin fluctuations. 

Our experiments were performed using triple-axis spectrometers at the High-Flux Iso- 
tope Reactor (HFIR) of the Oak Ridge National Laboratory [0,|J and the MARI di- 
rect geometry chopper spectrometer at the ISIS pulsed spallation neutron source of the 



Rutherford Appleton Laboratory (T^|. The momentum transfer (q x ,qy,qz) is measured in 
units of A -1 and reciprocal space positions are specified in reciprocal lattice units (rlu) 
(h, k, I) = (q x a/2n, q y b/27r, g 2 c/27r), where a, b (~ a), and c are the lattice parameters of the 
orthorhombic unit cell of (123)06+x- For this study, we prepared four single-crystal samples 
by the melt texture growth technique MM. Subsequent annealing resulted in oxygen stoi- 



chiometries of x = 0.6, 0.7, 0.8, 0.93 with superconducting transitions at T c = 62.7, 74, 82, 



92.5 K, respectively. 

The complete magnetic excitation spectra (Fig. 1) are obtained by integration of the 
magnetic neutron scattering over the two-dimensional Brillouin zone (BZ) of (123)06+a; for 



x = 0.6 at several temperatures Because the CuC>2 planes in (123)06+2; actually 

appear in coupled bilayers, magnetic fluctuations that are in-phase (acoustic or x" c ) or out-of 
phase (optical or x'L>) EH with respect to the neighboring plane will have different spectra 
EhE!- That the optical excitations occur at higher energies than the acoustic excitations, 
and are characterized by a substantial gap even in the normal state at 80 K, demonstrates 
that the bilayer coupling remains antiferromagnetic on going from the insulating parent 
to the superconducting metal. Upon lowering the temperature from 290 K to 35 K, the 
most dramatic change in the local-susceptibility x"{ u ) f° r energies less than 100 meV is 
the appearance of a resonance at 34 meV in the acoustic mode as shown in the shaded 
areas of Fig. 1. For x = 0.6, the resonance is clearly present not only below T c but also 
at 80 K, albeit in reduced form. Beyond the resonance, there is a broad feature extending 
at least to 220 meV, with a maximum at approximately 75 meV. This broad feature is 
reminiscent of the continuum, peaked at around 20 meV, seen in the single layer compound 



Lax 86 Sr .i4CuO 4 fjTBJ] . Well below the resonance peak, scattering is suppressed and a true 
spin-gap with a value of about 20 meV develops in the superconducting state (Fig. 1A). 
At room temperature, the optical and acoustic spectra (Figs. 1C and IF) show no maxima 
over the range of our measurements. 

Because the most dramatic change in the magnetic spectrum is the change in the reso- 
nance intensity, we focus on the temperature and composition dependence of the resonance. 
The detailed momentum and frequency dependence of S res (q,uj) obtained at HFIR at vari- 
ous temperatures around the resonance energy for (123)06.6 (Fig. 2) show that in the low 
temperature superconducting state (Fig. 2A), the spectrum is dominated by the resonance 
at 34 meV. At temperatures just above T c (Fig. 2B), the resonance broadens and decreases 
in intensity, consistent with the ISIS data of Fig. 1. The resonance peak intensity appears 
to shift smoothly toward higher energies in the normal state with increasing temperature. 



A signature of the resonance remains at 125 K and 150 K, but essentially vanishes at 200 
K. Previous temperature dependent measurements for underdoped (123)0 6+a; indicate that 
the resonance peak intensity changes its characteristics at T c @,[|]. While the intensity at 
the resonance energy has a noticeable upturn upon cooling through T c , the scattering at 
frequencies above and below the resonance decrease below T c ||||. This behavior can now 
be understood as arising from the narrowing of the resonance in energy from the normal to 
the superconducting states. 

The temperature dependence of the resonance peak intensity [S res (q,u)] for three 
(123)06+ie samples with different doping levels and transition temperatures is shown in 
Fig. 3, A throught C. For x = 0.6, the temperature dependence of the momentum- and 
frequency- integrated resonance is found to be different from that of the peak intensity (Fig. 
3A). For (123)06.8, similar behavior also occurs, but because the resonance is weaker above 
T c , the counting times required to obtain reliable integrals of the type found for x = 0.6 
are prohibitive, and so we are content with simply plotting the peak intensities. Finally, for 
x = 0.93, no broadening in either energy or wave vector has been identified and so in 

this case, the temperature dependence of the integrated spectral weight does actually follow 
the peak intensity. We define the mean-squared (fluctuating) moment associated with the 
resonance as {m 2 es ) = 3/(27r) / d(huj)Xre S (^) I '0- ~ exp(—huj/kT)), where % is the Planck's 
constant divided by 2n, k is the Boltzmann's constant, Xres( u ) 0jEOI is the resonance part 
of the acoustic spectrum of Fig. 1, and the factor of 1/2 arises from averaging sm 2 (q z d/2) 
over q z . For (123)0 66 and (123)0 6 . 93 our measurements are expressed in absolute units 
obtained by scaling to the low temperature measurements performed on both compounds 
at ISIS |TJ. 

The most obvious feature of Figs. 3, A throught C, is that as the doping level and T c 
decrease, there is a progressively larger pretransitional regime above T c . Specifically, for 
the ideally doped sample (x = 0.93), the onset of the resonance occurs at a temperature 
T* which almost coincides with T c . For x = 0.8 and 0.6, T* increases to approximately 
115 ± 15 and 150 ± 20 K, respectively, even while both T c and the resonance energy itself 
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are reduced. Thus, the weight of the temperature-dependent resonance joins the long list 
of properties which show pretransitional behavior in suboptimally doped (123)06+2. The 
cross-over temperature T* (Fig. 4A) coincides with the temperatures below which the 
temperature derivatives of the electrical resistivity dp(T)/dT |T9|J20|] and the Cu nuclear 



(T\T) 1 relaxation rate [21 reaches a broad maxima. The anomalies occurring at T* are 



generally associated with the opening of a pseudogap in the low-energy spin excitation 
spectrum, a supposition also supported by neutron scattering data such as those in Fig. 1. 

According to thermodynamics, a metal undergoes a transition into the superconducting 
state because such a transition can lower its total free energy, F [p2 |. The difference in free 
energy of the system between the normal state, extrapolated to zero temperature, (F N ) and 
the superconducting state (F s ) is the condensation energy, that is, E c = (F N — F S ) T=0 . In 
principle, the free energy of a system, and therefore Ec, can be derived from the temperature 
dependence of C el (x,T) ||. For (123)06+2;, C el (x,T) has been measured by Loram et al. 
||. The key features of the corresponding data (Fig. 3E) are a sharp jump at T c for 
the optimally doped (123)0 6 . 92 , which becomes much suppressed upon reduction of the 
oxygen content. Although there is less entropy released at T c , more seems to be released 
at temperatures above T c , as x is reduced below its optimal value. Thus, the specific heat 
tracks the temperature derivative (Fig. 3D) of the spectral weight of the resonance. Just 
as for the specific heat, decreasing doping reduces the maximum in d(m^. es )/dT at T c and 
introduces progressively broader high-temperature tails. 



In the t-J model |2j|, the Hamiltonian of the system consists of a nearest-neighbor hop- 



ping (t) term accounting for the kinetic energy (Eg) of the carriers and a magnetic exchange 
(J) term of the type describing the insulating antiferromagnetic parent compound ||24|| . Mag- 
netic fluctuations with a particular wavevector q contribute to the exchange energy (Ej) via 
a product with J(q), the Fourier transform of the exchange interactions [[K]]. Although 
magnetic spectra as complicated as those in Fig. 1 will be difficult to describe within such 
a simple model, the exchange part of the t-J Hamiltonian does provide a convenient way to 
quantitatively estimate the effect of magnetic fluctuations on the thermodynamic properties 



of high-T c superconductors. For (123)06+2, which has two coupled CuC>2 layers per unit 
cell, the exchange interactions consist of the nearest-neighbor spin-spin coupling J in the 
same CuC>2 plane and the coupling J± between two CuC>2 planes within the unit cell. Since 
J is much larger than J± , the exchange energy is (T^|11|J2^| : 

Ej * j^O 2 l^rl dq * dq *\ j[cos(g - a) + cos ^ b)] [Xac(g r-^p(-^/ ( y ' " )] ^ 

where g is the Lande factor (~ 2) and //£ is the Bohr magneton. Thus, Eq. 1 gives the 
contribution of the exchange energy to the total energy. Since the heat capacity is the 
temperature derivative of the total energy, we can estimate how the magnetic fluctuations 
responsible for the resonance affect the thermodynamic properties of (123)0 6+x ||||. As the 
resonance grows, more spins fluctuating at the resonance frequency are correlated antiferro- 
magnetically, and the exchange energy is correspondingly reduced. Assuming that the spins 
accounting for the resonance are not spatially correlated at high temperatures and therefore 
do not contribute to Eq. 1, the contribution of the exchange energy of the resonance {E r f s ) 
to the specific heat, or Cj es (x,T), is then 



fires/ ^ .V _ _ 7 / _V )_ AresV ^ I _ ^V__res£ fn\ 

J 1 ' J ~ dT 4fi 2 B dT\J n l-exp(-hu/kT)J 2/i| dT Kh 

Therefore, Cj es (x,T) is approximately proportional to the temperature derivative of the 
spectral weight of the resonance, as we already discovered from comparison of our data with 
the measured specific heat. In deriving Eq. 2 from Eq. 1 we relied on the fact that the 
resonance is sharply peaked around (n/a,n/b). The contribution of the resonance exchange 
energy to the heat capacity calculated using Eq. 2 with the measured exchange coupling 
J = 125 ± 20 meV for (123)06.6 13 is shown in Fig. 3D (right hand scale) and can be 
compared to the measured C el (x,T) of Loram et al. in Fig. 3E 0. The model calculation 
with no adjustable parameters shows that the resonance can provide enough temperature- 



dependent exchange energy to yield the superconducting anomaly in the specific heat |26 
This suggests that a large part of the specific heat near T c and its associated entropy are 
due to spin fluctuations. 



Clearly, the measured specific heat includes contributions other than that due to the 
exchange energy of the resonance (Ej es ). While we included the most obvious temperature 
dependent feature in our calculation, the temperature evolution of the remainder of the 
spin excitation spectrum is neglected. Most notably we ignored the formation of the spin 
pseudogap below T* and spin gap below T c , both of which could increase Ej. Further, 
we have not considered the contribution of Ek to C el (x,T). For superconductors, Ek is 
expected to increase from the normal to the superconducting state JTOf. Both this effect 
and spin (pseudo)gap formation could partially compensate for the reduction of Ej due 
to the appearance of the resonance. Such compensation might explain why the area of the 
calculated Cj es (x, T) anomaly is considerably larger than that of the C el (x, T) measurement. 

In the next few years, the contributions missing from the present analysis should be- 
come calculable using better neutron and optical conductivity data, which would provide 
the necessary information on the temperature-dependent magnetic exchange and kinetic en- 
ergy terms, respectively. It would also be of considerable interest to make further model 
calculations to estimate the change in the resonance exchange energy between the normal 
and superconducting states at zero temperature and hence the resonance exchange energy 
contribution to the condensation energy ]IUyTT |. Unfortunately, this comparison is difficult 
because the normal state S(q, uj) must be extrapolated to zero temperature. This cannot be 
done reliably with the present data and presents a challenge for the future. 

In summary, we show that a pseudogap regime bounded by a temperature T* > T c deter- 
mined from transport, NMR, and thermodynamic measurements can be associated with the 
appearance of the resonance peak and the suppression of the low frequency response in the 
magnetic excitation spectrum. Our complete spectra allow us to establish the weight of the 
resonance relative to that of other spectral features. A simple calculation of the exchange 
energy using the measured temperature dependence of the resonance shows that spin fluctu- 
ations can account for a large part of the electronic specific heat near the superconducting 
transition. 
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FIGURES 

FIG. 1. Local- or wavevector- integrated frequency-dependent magnetic susceptibility x"( w ) for 
acoustic and optical modes in (123)0"6.6 corrected for the isotropic Cu 2+ magnetic form factor and 
instrumental resolution [13]. The effect of Cu 2+ anisotropic form factor [S. Shamoto, M. Sato, J. M. 
Tranquada, B. J. Sternlieb, and G. Shirane, Phys. Rev. B 48, 13817 (1993)] is not considered. The 
data were collected with the crystal mounted in the (h, h, I) scattering plane at MARI spectrometer 
at ISIS. (A-C) were chosen at q z positions to emphasize acoustic modes and (D-F) the optical 
modes. Integrating the area under the resonance peak and converting to moment squared yields 
(m 2 res ) = 0.06 ±0.04 //|f.u. _1 at 35 K and 0.04 ±0.03 ^f.u.^ 1 at 80 K. The resonance in (123)0 6 . 6 
accounts for about 14% of the spectral weight up to 80 meV at 35 K (Fig. 1A) and about 1% 
of the total moment-squared 2g 2 y? B S(S + 1) = 6 n 2 B expected per formula unit (f.u.) assuming 
S = 1/2 for each copper ion. Similarly, we find that the 41 meV resonance for (123)06.93 has 
( m res) = 0.06 ±0.04 /Xgf.u -1 at 25 K. Note that the estimated error in (m 2 es ) does not include the 
~30% systematic error involved in normalizing the scattering to a vanadium standard. 

FIG. 2. The magnetic structure factor S res (q,oo) around the resonance energy for (123)06.6 at 
various temperatures. The data were taken with the HB-3 triple-axis spectrometer at HFIR using 
pyrolytic graphite as monochromator, analyzer, and filters. The collimations were, proceeding 
from the reactor to the detector, 50'-40'-80 / -240 / , and the final neutron energy was fixed at 30.5 
meV. To compose the color figures, constant-energy scans along the (h, 3h, 1.7) direction [3,4] 
from h = 0.42 to 0.58 rlu were performed at energy transfers of fiuj = 25, 27, 30, 34, 37, 40, 
43, and 46 meV. Magnetic intensity was assumed to be scattering above the linear background 
of the constant-energy scans. The scattering at different temperatures was normalized to 300 
monitor counts which correspond to 7 minutes per point at Twj = 34 meV. Note the change in the 
vertical color bar scales, made to better display the dramatic decrease of the magnetic intensities 
on warming. 
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FIG. 3. Temperature dependence of the resonance peak intensity (filled circles) for (A) 
(123)0 6 . 6 at hu = 34 meV, (B) (123)0 6 . 8 at 39 meV, and (C) (123)0 6 . 93 at 40 meV at the 
(0.5,1.5,1.7) rlu position on triple-axis spectrometers at HFIR. Left axes show intensity normalized 
to (rn* es ) and right axes to the peak intensity at low temperatures. The error bars in (A-C) 
include only the statistical errors of the triple- axis measurements. The squares in (A) are the 
integrated S res (q,u>) and the triangles in (C) are from Ref. 2. T*s are defined as the initial ap- 
pearance of the resonance. The solid line in (C) below T c is a fit using the modified two fluid model 
[/(T)/J(0) = 1-(T/T C ) 3 - 09 ]. (D) The estimated d(m 2 res )/dT for (123)0 6+a; is from the temperature 
derivative of the solid lines in (A-C). To calculate the absolute magnitude of Cj es (x,T) using Eq. 
2, we have assumed J = 125 meV (Ref. 13), (m 2 es ) = 0.06 ^|f.u _1 , and 1 mole = 666.15 grams for 
each of the three compositions. Note we have not considered the uncertainties of Cj es (x, T) arising 
from the systematic errors in determining the absolute magnitude of (m^ es ), due most notably to 
the ambiguity in removing the nonresonance portion of the signal (Fig. 1A). (E) C el (x,T) from 
Ref. 9 was converted to SI units using 1 g-at.= 1/(12 + x) moles. 

FIG. 4. Phase diagram of (123)06+x summarizing the results of transport, NMR, and neutron 
scattering. Because various methods of oxygenation, slightly different T c values are reported for 
the same nominal doping of (123)06+a;. We avoid this inconsistency by plotting characteristic 
temperatures as a function of T c . (A) The open circles and open squares are temperatures which 
dp(T)/dT reaches a broad maximum (Refs. 19, 20). The open diamonds show the the pseudogap 
temperature T* determined from NMR measurements (Ref. 21). The filled circles locate T c and 
T*, where the resonance first appears in our samples. (B) Filled circles show resonance energy 
versus the transition temperature T c . Filled squares are from Fong et al. (Ref. 5). Horizontal 
error bars are superconducting transition widths. The solid lines are guides to the eye. 
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